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TG: Distriljution List 

FEOMj J. Uskavitch 

SUBJECT; MIG Subroutine Library 

Group 22 Internal Memorandm Ho. 22“SA-56 


Abstracts FiTeo general-utility subroutines are presented in coded form 
for use on MTC. Information is also given about average oper¬ 
ation timOj, accuracy, storage requirements and operating de¬ 
tails pertinent to each subroutine. Explanatory notes are 
also included with each one for the aid of programmers. 

The subroutines included are; 

^ a^^ Addxtxon, Wxth Scale Factor .o.o.oooooo.ooo .Page 3 
^b) [Decimal Prxntout o.......©..'...............©Page ^ 

(c} Point Bxsplay ooooooooooooo.©.............©Page 11 

(d) Axes—Grid Bxsplay .©.©©©©o©©©. ©©©©©©.©©©©©©©Page 

(e) Special Square Boot©o©©©©©.©©.©©©©©©.©©©•©.©Page 19 

A diagram at the end of the memorandnm shows one set of storage 
registers for eight subroutines for which tapes are available© 


General Remarks; 

All of the subroutines given in this memoraadtim have been test¬ 
ed on MTG and are believed to be free of errors© However, if. mistakes 
are found in any eases, the author will sincerely appreciate having his 
attention called to them© 

In developing these subroutines, the writer has attei^ted to 
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make them as easy to use as possible. Consequently^ time and space have 
sometimes been sacrificed to increase the ease of programmingj moreover, 
whenever possible, space has been saved in preference to a few miero-= 
seconds of operating time. However, the writer will sincerely appreciate 
receiving any information concerning excessive use of either space or 
time in these subroutines. 

As will be noted in the following pages, all addresses have 
been written as octal numbers, thus minimizing the work needed for trans¬ 
cription into the normally used tapes. However, in all other respects, 
the notations used are identical with those of Memorandum 6M“3^9r 
eluding the re-use of the letters a,b,e,d, and e to denote the individual 
subroutines. The waiter has reused these denotations because he did not 
know what other additions were being made to the subroutine library and 
consequently did not know just where these subroutines would fit in the 
library. In any ease, when the aimber of MTC subroutines become large 
enough, he hopes that Group 62 will re-issue all of them with identical 
formats in one memorandum. 

Tl® duration time--givaa for each subroutine must be considered 
as only a rough approximation. All operations have been assumed to 
occur within core memory, and individual execution times were obtained 
from Drawing B-62366, which follows page 30 of 6M-2527-2. However, in 
all of the subroutines, the exact operating time depends on the input, 
and only an approximate averaging was performed overVthe possible inputs. 

The accuracies associated with these subroutines are all de¬ 
rived from round-off errors and were calculated in the same manner as 
described in 6M“3^97 « ®3-e one exception to this statement is found in 

subroutine (e), the accuracy of that output being determined in part 
by that of SRL #r of 6 m-31i-97. 

The general notation used in this memoranda may be summarized 
as followss 

pN “ relative address of a register in the operating part of 
a program, where p is a lower-case letter designating 
a particular subroutine and H is an octal integer. 

Kpl g relative address of a register containing a constant 

used in the program or being used for temporary storage 
where p and N are as above. 

k s the field in which the particular subroutine operates, 

m s the field in which the main program operates. 

Included at the end of this memorandum is a diagram which de¬ 
picts the storage locations utilized by available tapes of several of 
the MTC subroutines. The writer gladly offers to lead either the flexo 
or the converted tape if a reproduction is desired by anyone. 
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(a) ADDITION, WITH SCALE FACTOR 


1 1 

Input : Store the nuaiber (N^ 2 ) of the numhers (x;^)to he added in 

Ka3 and the extended address of the first number (x^^) in Ka4. 
Enter at aO after ra al30 with m. Content of AC is imiaater- 
ialo (All Xj^ must have been previously stored in conseeutive 
registers in a single field.) 


Output; 


2-n ^ 


x^ (the 15 most significant binaary digits of 



in AC# with the scale factor n left in Ka5. (The remaining n 
digits of the sum will be left in the first n digit positions 
of BRj i.e., positions 0 through n-1.) If the result is zero, 
it will be a negative zero, unless all x a +0. 


Other Subroutines Needed; None. 


Duration ; Approximately 915(N4'2) microseconds. 

Jl . 

Accuracy; Maximum error in ^ *i 2 for a>0, (For n s o, no error 
exists.) The answer obtained is exactly that which one would 
gat by just taking the 15 most significant binary digits of the 
stm without first roxinding off. Therefore, the results in the 
AC is always smaller than the true value, within this toler¬ 
ance. No error, if the n digits in the BR are also utilized. 


Storage; 90 program registers plus 9 constants or temporary storage. 

Notes; (l) The original numbers (x^^) are not destroyed by the addition 

process. 

(2) If x^ is in register 37^ of field 3# the extenied: address 
should be stored as O. 30376 . 

( 3 ) The scale factor used is always the smallest one possible. 
For example, as long as the sum requires no more than 

15 digits, no scale factor is used (n « O). 
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(a) 


ABDITIOU , WITH SCALE FACTOR 


Programs 


(enter) aO 

rf 

al31 

set up to leave subroutine 

1 

ca 

KaO 

+0 

. 2 

st 

Ka5 

a 

3 

St 

Ka6 

2; n 

h 

St 

Ka7 


5 

ca 

allj- 

srlk 

(a53)^ 6 

ra 

ali.5 


T 

ca 

al4 

srl4 prepare for 

alO 

su 

aii-5 


11 

ra 

a23 

- J 

12 

ca 

Kah 

extended address of x. 

13 

ra 

a21 

Ga(RCx2^) 

14 

sr 

14 

15 . 

ra 

a20 

sof( ) 

l6 

cs 

Ka3 

K 

(a37)-» 17 

sis 

KalO 

counter 




initial 


conditions. 


a20 

sof— 


21 

ca — 

^i 

22 

sof 

k 

23 

sr ■ 

— 

(0,3,6,11,: 

2k 

et 

Ka2 

1,00007 

^ 25 
: 26. 

tn 

a27 


^r 

a30 


27' 

(su 

Ka2 

1,00007 

-a30 

^ad 

Ka6 

^ n 

31 

St 

Ka6 

^ n 

32 

ca 

a21 


33 

ad 

Kal 

41 

3^ 

ra 

a21 


35 

ca 

Kal 

+1 

36 

sm 

KalO 


37 

tn 

al7 

have these 
if no, reti 

aUO 

ca 

Ka6 

Zn 

41 

sr 

3 


42 

st 

Ka6 

Zn 


obtain x. in AC. 

1 

extract tbe sign plus three 
/ digits, beginning with the 
J 3 most insignificant 

is the sign negative? 
if no, transfer to add. if 
yes, correct extracted'number 
to provide true value of di- 
X ts o 

1 add to sum of same digits 
T obtained from preceding x's 
3 and stote. 

") prepare to extract 

1 ^ 

add-f 1 to counter. 


if yes, shift 3 most insigni- 
fieent digits of the s\am into 
the BE. Store those left, 
for addition of next 3 ex¬ 
tracted digits. 





22-SA-56 


Page 5 


Programs 

(a 

) ADDITION, 

a43 

sr 

15 

44 

er 

403 

45 

sr 


46 

ad 

KaT 

47 

St 

KaT 

a50 

ea 

a45 

51 

su 

a4l 

^52 

/ 

tn 

a54 

. / 53 

tr 

a6 


WITH SCALE FACTOR 


55 

56 


ca KaO 
am Ka6 

tn a6l 


( 14 , 11 , 6 , 3 , 0 ) 

rr i 

Hv 


obtain ttie remainder in the 
AC in correct position for 
addition to preceding remain¬ 
ders » 

add to preceding reminders 
and store. 

sr( 14,11,6,3,Op prepare to change the posit-, 
sr 3 |lon in the AC of the next 

jremainder. 

have all digits "been added? if yes, 
to combine the tvo stuns £.nrand;^r. 
return to reset initial conditions, 
add next 3 digits. 

if yes, 15 most insignificant di¬ 


transfer 
If no, 
and to 


-hO 

Zn 


1 57 

ca 

Ka7 

ITr 

^ a60 

tr 

al30 


^X6l 

ca 

KaO 

40 

62 

sm 

Ka7 

2Lr 

/63 

tn 

367' 


/ 64 

1 65 

ca 

Ka6 


sr 

21 


1 66 

St 

Ka7 


67 

es 

Ka6 

r n 

a70 

lah 

Ka7 

r r 

71- 

tn 

al06 



72 

ca 

Ka6 

Z n 

/73 , 

tn 

a76 

/ 74 

^ v su 

Kal 

4-1 

V 75 

>tr 

a77 


476 

( ad 

Kal 

A-1 

77 

Vst 

Ka6 

Z n 

alOO 

ca 

Ka7 

X r 

,^101 

tn 

al04 


/ 102 

su 

Kal 


1 103 

✓tr 

al05 

^104 

( ad 

al 

4- 1 

105 

>st 

TSa.7 

£ ^ 


1 





gits are now inEr, with the larg¬ 
er ones being in ^ n. 
is|^ nj - 0 ? 

if no, transfer to continue* 

if yea, putSr in the AC, and 
transfer out. 

ifl^nj ^0, 
isl Z r I 0? 

if no, transfer to continue. 

if yes, ma^e^r the same sign as 
^ n. (this particular amount of 
of shifting was picked to also 
provide a needed constant.) 


are the signs of^n and^r the same 
if yes, transfer to find n, and to 
combine, 

if no, proceed to prepare the two 
sums for proper eambiningb 
decrease the magnitude 0fS"n by 1, 


decrease the magnitude ofZr by 1. 
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(a) ..ADDITION, WITH SCALE FACTOB 


Program: 


(a7l)->106 

ca 

Ka6 

107 

9M 

BsQ) 

( -jiallO 

Bt 

KalO 

/ 111 

8U 

Kal 

I 112 

tn 

al21 

113 

ca 

Ka5 

114 

ad 

Kal 

\ 115 

8t 

Ka5 

\ 116 

ca 

KalO 

\ 117 

sr 

1 

'al20 

tr 

alio 


S. a 
-^0 
|Snt 

+ 1 


to find n, first ottain |£ln| . 

iajsnl , or subsequent remainder 
= 0? if yes, transfer to combine to 
Zn and^r. 



if no, add 41 to the scale factor 
n. 

shift 1 digit out of jS'nl 


transfer back to store the remainder and to 
ask if it is zero. 


(all2)-4121 

ca 

Ka5 

122 

ad 

065 

123 

ra 

al27 

124 

ca 

Ka7 

125 

sr 

17 

126 

ca 

Ka6 

127 

cr 

— 

(a60)»al30 

sof 

mi «■» 

131 

tro 



n 

sr21 


prepare to obtain the 15 most sig¬ 
nificant digits in the AC. 


] 


(214N) 


put.^r in digits O-I5 of BK. 

put n in the AC. 

obtain the I5 ^st significant 

digits in AC. 




leave subroutine 


e 


KaO 

+0 



1 

-hi 

.*.. " ' 


2 

1.00007 

extractor 


3 

— 

N = no. of numbers 

to be added 

4 

5 

—- 

extended address of x. 
a = scale factor 

6 

— 

2:n 

1 

7 

— 

^ 11 

f temporary storage 

KalO 

— 

counter; |^a| -40. 

) 
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(b) DECIMAL PRINTOUT 

Input: Store the position (n) of the binary point in Kb30: n = +0 

for an integer, +0<n<4- 15 for a non-integer, n - +15 for 
a fraction less than onCo After ra b53 with m, enter with 
X in AC either at bO (if sign is wanted) or at blO (if only the 
magnitude is desired)o 

Output: Prints + or - sign, if wantedo Then, if x is an integer, 

prints five decimal places, with initial zeros being replaced 
by spaces; if x is a fraction, then prints a decimal point 
followed by four decimal digits; if x is a non-integer >1, 
prints the integral part followed immediately by the decimal 
point and fractional part (each as described in preceding 
cases)o Two spaces are printed after every numbero 

Other Subroutines Needed: NonCo 

Duration: approximately 0ol25 second for the sign plus an average of 
Oo 9 seCo for an integer or fraction less than one, and lo 8 
seconds for a non-integer greater than onCo 

Accuracy: No error for an integero A maximum error of Oo 0001 

(decimal) for both of the other caseSo The fraction printed 
is exactly that which one would obtain by taking only four 
decimal places of the true value without rounding off; there¬ 
fore the fraction printed is always lessj than its true value, 
within this toleraneeo 

Storage: 77 (decimal) program registers plus 27 constants or temp¬ 

orary storagCo 

Notes : If the integral part of x is known to have less than five sign¬ 

ificant figure s, time and space can be saved by ra Kb 24: with 
(or st Kb24^) the address of the register containing +1000, 
+100, +10, or +1 depending on the number of significant fig¬ 
ures* 

2, For those more used to the conventional scale factor, it may 
be well to note that the position of the binary point is ess¬ 
entially the same as (+15- n), when the’’scale factor is 
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(b) DECIMAL PRINTOUT 


Program 
(enter if- 
sign is 
wanted.) 


bO 

1 

' 3 
4 
>5 
6 
7 


b54 

Kb31 

b5 

Kbl7 

b6 

Kb20 

400 

bl2- 


X 


-1-13(-!-) 

+29(-) 

J 


set up to leave subroutine 0 


) find sign of x and print- 


(enter if—>bI0 

f 

b54 


only mago 11 

1 St 

Kb31 

X 

is wanted) 12 

^ cs 

Kb 31 

X 

/13 

tn 

bis 


' 14 

St 

Kb31 


ills 

cs 

Kb30 

n 

16 

ad 

Kb 7 

+1S 

17 

tn 

b37 


bZO 

su 

Kb 7 

+1S 

21 

^-tn 

/ 

b24 


22 

1 

, tr 

1 

bS5 


(bl02) —23 

! tr 

bSO 



(blOZ) 
{bl7) . 


24 

^ cr 

2000. 


2S 

ca 

Kb30" 

n 

26 

ra 

b30 

sr(n) / 

27 

ca 

Kb31 

|x| 

b30 

sr 

(n) 


31 

St 

Kb31 

x.-^x| j 

32 

cr 

420 

1 

33 

St 

Kb32 

0 X ’ r 

34 

tr 

bSS 

_ ) 

> 3S 

ca 

Kb32 

. X ”*7 

36 

st 

Kb31 

. x-^xjj 

- 37 

ca 

Kb21 

+ 17 ( 0)7 

b40 

pr 

400 

J 

41 

ca 

Kb 31 

/xi 7 

42 

rnh 

KbO 


43 

st 

Kb31 

|x» -» 

44 

ca 

Kb23 

suKbB 

4S 

ra 

b6S 

T 


obtain and s|:ore magnitude 

of 


does n = 15? 

if yeSj, transfer to print a 
fractiono 

if no, does n = 0? 

if n^l 0, transfer to print 

non-integer^ 

if n=0, transfer to print an 
integer» 

return, and transfer to leave 
sub routine p 
clear BR 

obtain and store for printing 
the integral part of x. 


obtain and store temporarily 
the fractional part of Xp 
transfer to print an integer 

return, and prepare to print 
fraction. 

print decimal pointo 

find 10, 000 times the fraction 
and prepare to print as integer 


pare to start with a 
subtraction of +1000. 
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(b) DECIMAL PRINTOUT 


46 

ca 

bill 

tn b76 

set up to print initial zeros. 

47 

tr 

b60 


transfer to print. 

(b23,bl02)-»b50 

ca 

Kb22 

+8(spac^ 


51 

pr 

400 

1 

(when X has been printed. 

52 

pr 

400 

J 

'print two spaces. 

53 

sof 


1 


54 

tro 

-- 

i 

leave subroutine 

(b22,b34)—»55 

ca 

Kb24 

su KbO') 

when an integer, prepare to 

56 

ra 

b65 

J 

start with subtraction of +10000 

57 

ca 

Kb25 

tn bl07? 

prepare to replace initial zeros 

(b47)-^ b60 

ra 

b66 

J 

with spaces. 

'61 

rf 

bl02 


set up to return. 

(bl06)-^ 62 

ca 

Kb26 

ca Kb57 

initially set flexo code 

, 63 

ra 

b76 

J 

to print a zero. 

(b75) -^ 64 

ca 

Kb31 


is |x| or later remainder^ 

65 

su 


(Kb0-Kb4) 

10, 000 (1000,100,10 or 1)? 

66 

tn 

-- 

(b76 orblO?') if no, transfer to print the 





digit or space if initial zero 





not wanted. 

67 

st 

Kb31 

jx 1 

if yes, store the remainder. 


b70 

71 

72 

73 

74 

75 


ca 

ad 

ra 

ca 

ra 

tr 


b76 

Kb4 

b76 

bill 

b66 

b64 


(b66,blll) 


(bll4) 


(b66) 


► 76 

ca 


77 

pr 

400 

blOO 

ca 

Kb27 

101 

su 

b65 

102 

tn 

-- ( 

103 

ca 

b65 

104 

ad 

Kb4 

105 

ra 

b65 

106 

tr 

b62 


+1 f add 1 to the Flexo-code 
storage, to be printedo 
tn b76 1 since this digit is^O, prepare 
J to print subsequent zeros, 
_ return to subtract again, 

:Cb5-Kbl6)y obtain Flexo-code in AC 
Jand print the digit 


su Kb3 
su Kb0-Kb4); 


is this the last digit? 
if yes, leave. 


7 


+1 


if no, prepare for subtraction 
of next smaller power of 10, 
return to find next digit# 


^'.107 

ca 

Kb27 

bllO 

su 

b65 

111 

tn 

b76 

112 

ca 

Kb22 

113 

pr 

400 

114 

tr 

bl03 


suKbS /if digit is an initial zero, 
suKb0-Ko4 is it the last digit? 

—^if yes, transfer to print it, 
+8(spac^ if no, print a space, 

_prepare to find next digit. 
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(b) DECIMAL PRINTOUT 


KbO 

+10000 


1 

+1000 


2 

+100 


3 

+10 


4 

+1 


.5 

+62 

(0) 

'6 

+21 

(1) 

7 

+15 

(2) / 

KbIO 

+7 

(3). / 

II 

+11 

(4) 

12 

+19 

(5) V 

13 

+27 

(6) . 

14 

+23 

(7) / 

15 

+3 

(8) 

16 

+54 

:{9) \ 

17 

+13 

(+) _J 

Kb20 

+29 

(-) 

21 

+17 

(.) 1 

22 

+8 

(space)f 

23 

su Kbl 


24 

su KbO 


25 

tn bl07 


26 

ca Kb 5 


27 

su Kb 3 


Kb30 

-- 

n { postic 

31 

-- 

X, jKl , 

32 


, x{fracti( 


Flexo codes 


Flexo codes 
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(c) POINT DISPLAY 

Input ; Having ra c36 with. m,put abscissa (x) in KclO and ordinate (y) 
in Kell and enter at cO. Content of AC is imnaterial. (The 
constants required in registers Kcl - Kc7 and the two er in¬ 
structions in c21 and c25 must have previously been stored. 

See notes below.) 

Output ; If —^max Fmln -y—ymax -^ where the minimum 

and maximum values are the desired limits for plotting, the 
point (x,y) will be displayed N times on the oscilloscope with 
whatever linear scale desired. If either coordinate is out¬ 
side the set limits, no point is displayed and the program 
may be halted. (See notes below.) 

Other Subroutines Needed ; Nona. 

Duration; (413 ■/* 76 n) microseconds. 

Accuracy: Each coordinate is displayed with an accuracy of 1 part in 

2 ^. 

Storage ; 35 (decimal) program registers, 11 (decimal) constants or 

temporary storage. 

Notes; 1. Since the intensity of the displayed point may want to be 

varied—either to distinguish one set of points from another 
or to provide for easy distinction from a superimposed grid, 
the nxunber of times (N) that the point is to be displayed must 
be stored in Kel, 

2, When the desired coordinate limits have been determined, the 
coordinates of the middle point (xm,yjjj) must be stored in Kc2 
and Kc5> respectively. If this point is displayed, it appears 
at the center (40,'^0) of the display. 


Examples: 

^ax 

■+.9999 

+200 

+20000 


■%in 

-»9999 

+100 

-20000 

2(W^ymin )• 

Ke2 

40 

+150 

+ 0 


Y 

“'max 

4100 

— 1 

+-«9999 


^min 

•40 

-11 

+.5000 

si yka^c+ymin) s 

Kc5 

450 

-6 

4,7499 


3. Since the limits of the display may be smaller than those of 
the available data, the desired ranges of the coordinates 
(with a factor of ■§•) must be stored in Kc3 and Kc6, The 
factor of y is included to provide for the times when the range 
is greater than 0.77777■> A continuation of previous examples 
gives: 
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(c) POINT DISPLAY 

i(%ax ” %Ln); Ke3 -^9999 +50 4-20000 

l(ymax - Kc6 +50 +5 +.2499 

4. To provide for a display with the proper scales the two scale 
factors must he stored in Kc4 and Kg 7. To he safe^ each 
scale factor must he such that the product of it and the half- 
range^ +511| i.e., (Kc3) (Kc 4) <511 and (Kc 6 ) (KcT)^ 511. 
continuing the previous examples: 

x-seale factor; Kc4 +511 +10 +.0255 

y-scale factor; Ke7 +10 +102 +2044 

5. The er instructions to he stored in c21 and g25 are determined 
hy the ranges and scale factors used. If hoth are^+1 and 
short-cyele is used, 12 (octal) places must he cycledj if 
either is^+1, 31 (octal) places must he cycled. Concluding 
the previous examples; 

program x; e21 cr431 cr4l2 cr431 

program y: c25 cr4l2 Gr4l2 cr431 

6. If the point to he displayed lies outside the predetermined 
region, one of several actions may he wanted. As programmed 
here, ha 0 is given and pushing the restart button changes 
frames and again halts the program. However, if the point 
should merely he ignored, then the instruction tr 036 can he 
put in c4lj if one wished to return to a different part of the 
main program when this happens, an sof and a tro instruction 
can he stored in e4l and e42 respectively. 
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(c) POINT DISPLAY 


Program; 

(enter) eO rf c37 

1 ca KclO 

2 su Kg2 

3 to e40 
k St KclO 

5 ca Kc3 

6 sm KclO 

7 tn c40 

clO ea Kell 

11 su Kc5 

12 to c40 

13 St Kell 

14 ca Ke6 

15 sm Kell 

16 tn e40 

17 ca KclO 

g 20 mh Ke4 
21 cr ( ) 

__22__st__KclO 

23 ca Kell 

24 mh Ke7 

25 cr ( ) 

26 St Kell 

27 es Kcl 
•:ijc30 st Kel2 

31 ea Kell 
I 32 ds KclO 
33 ea KeO 
^ 34 sm Kcl2 
V 35 tn e30 


36 sof 

37 tro 


set up to leave subroutine. 


X 

^ " 2(-^max'*'^min) 

x(- x-xm) 

x/2=4(x -V . /) 
' 2\ -^irr 

I I 




obtain the difference be¬ 
tween X and the median of 
the range, if this differ¬ 
ence produces an overflow, 

X lies outside desired range^* 
transfer to halt. 

^ax * 

if not, transfer to halt. 


y 

yin 


1 


y(-y-yj 

Ay/2 

1 y-yj, I 





obtain the difference be¬ 
tween y and the median of 
the range, if this differ-=r 
ence produces an overflow, 
y lies outside desired 
ranges transfer to halt. 


is iy^ 7 

if not, transfer to halt. 



x-seale factor 


scale the distance from 
display center in x-direction » 


put in Alin position to display.(note5.) 


X“Xi 


•m 


y^ym 

y-seale factor 


scale the distance from 
display center in y-dir- 
eetion. 


put in AG in position to display. (Note 5o) 


- 1 : 2 a— 

N 

counter 

yym 

x-Xjn 

- 1-1 

counter 



set up counter, 
display (x,y)* 
add-{-l to counter. 


if not displayed N times, go back to display 
again. 


leave subroutine. 
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(e) 

POIHT DISPLAY 

chO 

(ha 0) 

if X or y is not in desired range, halt. 

kl 

(op 0) 

'l upon RESTART, change frame and halt. 

h2 

(ha 0) 

\ (Note 6.) 


KcO +1 
1 
2 
3 
I). 

5 

6 
7 

KelO 

11 

12 


H: no. of times each point is to he displayed 
^ = i(%ax+'%in)* ot desired range. 

Ax;/2 s (Xfflax thalf the desired range. 

X“Seale factor 

yjQ s i(ymax + ymin)* “®^ian of desired range. 

Ay/2 s ■|( ymav " ymin^* desired range, 

y-scale factor 
xj x-x^; ■-) 

. r temporary storage 
j} J 

counter 
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(d) MBS GRID DISPLAY 


Input t First ra 036 with m. Then enter at dO if a new frame is wanted 
in the camera, otherwise at dl. Content of AC is immaterial. 
(The constants required in Kdl-KdT must have "been previously 
stored. See notes helow.) 

Output; A set of orthogonal axes displayed M times on the oscilloscope 
with any predetermined origin. A set of grid lines, with any 
predetermined spaeings, displayed one® to provide linear scales 
measiired from the origin. (See notes below.) 

Other Subroutines Heeded; (e) Point Display. 

Duration; . To display axes; 

£ ( 13to)(no. of points / axis) + (13^) (no.of points/ grid linej 
microseco Maximum — 13T(JJ-*l) milliseconds. 

For each grid line displayed, add; 

Ql09-f 67 (no. of points/grid 11miJ microseconds. 

Maximum—68 milliseconds / line 

Accuracy; Ho inaccuracy. 

Storage; 63 (decimal) program register, 8 constants, and 6 registers of 
subroutine (e). 

Notes; 1. Although this subroutine req;uires only 6 additional reg¬ 

isters to make it independent of the preceding routine, 
the 'writer cannot,foresee an independent use and consequen¬ 
tly has saved 6 registers. The changes required to make 
this an independent subroutine are believed to be obvious. 

2. This subroutine may be much longer than actually necessary 
becaixse of the many variations permitted; not only May the. no. 
of times (h) that the axes are displayed be changed indep¬ 
endently of the no. of times that points are displayed in 
subroutine (e), but also the densities of the points in the 
axes and in the grid lines may be changed independently. 

This great variability has been included since large 
variations•in the display equipment and/or photographic 
equipment and developing have in the past made it impossible 
to choose optimxmi values. 

3. The densities of the axes and grid line® are determined by 
the values of ^^andAgput in K<32 and Kd3o These numbers 
are easiest stored as positive octal, numbers, the first 
three digits giving the spacing between points and the last 
two being zeros. (See Note 6 for an example.) 
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k. The coordinates of the origin are also most easily stored 
as octal numbers, remembering that the ^ instruction only 
uses the sign plus first three digits. If a coordinate 
is positive, the last two digits should be zeros, as is the 
case for the densities; if a coordinate is negative the 
last two octal digits shoxild be made 77- (See Note 6 for 
example.) 

5. Instead of storing the distance between adjacent grid lines, 
the half-distances are to be put in Kd6 and Kd7. Kd6 pro¬ 
vides the spacing of the vertical lines, while Kd7 gives 
that of the horizontal ones. The half-distance is used to 
provide for the case when the desired spacing is greater 
than 0.77700* Each half-distance can be stored as a pos¬ 
itive octal constant with a zero in the last digit, the first 
four digits producing the spacing. (See Note 6 for an 
example.) In case no grid lines are wanted, merely store 
0.777^0 in both registers. However, the axes are always 
further intensified by one pair of superimposed grid lines, 

6. Example; Origin to be located in center of left half-plane. 
The axes are to be displayed 5 times with only 50^ of the 
points being displayed. Only one out of every eight points 
is to be displayed in each grid line. The vertical grid 
lines are to be 100 points apart, whereas the horizontal 
spacing is to be I75 points. 

Kdl +5' Kd5 0.00000 

2 0,00200 6 0.06200 

3 0.01000 7 0.12740 

4 1.37777 (or 1 . 40077 ) 
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(d) 

Program; 

(eater,if 
frame to "be 


changed) 

dO 

op 

0 

if not, 

1 

rf 

c37 

enter) 

2 

cs 

Kdl 

(d25) — 

3 

St 

Kcl2 

4 

ca 

Kd4 

5 

St 

KclO 

6 

ea 

KdO 

71 7 

ds 

KelO 


/dlO 

ad 

Kd2 


i 11 

/to 

dl3 

\ 12 

' tr 
\ 

d7 

13 

ca 

KdO 

yflk 

St 

KelO 

/ 15 

ea 

Kd5 


16 

ds 

KelO 


17 

ca 

KelO 


1 d20 

ad 

Kd2 

\ 21 

/to 

d23 

^22 

' tr 

dl4 

23 

\ 

•^ca 

KcO 

2h 

sm 

Kcl2 

25 

tn 

d3 



26 

ca 

Kd4 

(<134)- 

27 

St 

KelO 


d30 

su 

Kd6 


/31 

to 

d35 


1 32 

i ^33 

su 

Kd6 


to 

d35 


\l 34 

am .V— — — 

tr 

d27 

(d50)- 

-^^35 

ca 

KdO 


-36 

ds 

KelO 


r 37 

ad 

Kd3 


1 d40 

to 

d42 


\4l 

tr 

d36 


AXES - GRID DISPLAY 


change frame. 

set up to return to main program. 

® 7 set up counter# 

coxmter J 


X 

1.00077 

X 


Aa j producing 


display 

y-axis 


1.00077 
X->0.77700 

yo 

X-^. 77700 
x-> 0.77700 


display 

x-axis 


[-add-t-l to counter, 
counter J 

if not displayed H bim^; go hack to display 
again e _ 




x-^X“positi6n of. 
first grid line. 
Ax/2 

Ax/2 


J 


find coordinate of 
first vertical grid line, 
position has been 
found -when overflow’ 
occurs. 


1.00077 

'^x “ coordinate 
[of grid line 

producing 






display Vertical 
grid line. 
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(d) AXES - GRID DISPLAY 


(dlK))" 


-^42 

ea 

KelO 

X 

\ find coordinate of next 

k3 

ad 

Kd6 

A x/2 1 

' vertical grid line. 

/ kk 

to 

d51 

\ 

1 all vertical lines 

/ 45 

ad 

Kd6 

^yil2 

/ have been displayed 

1 A6 

to 

d 51 


when overflow occurs. 

1 / 

V d 50 

i_ 

st 

KclO 

x4 A. X 


tr 

d35 

if vertical grid incomplete, return to dis 
play next line. 


^ 51 ca 
52 St 

f 

' 53 su 
5 ^ to 


d 70 ^c& 

71 ad 

72 to 

73 ad 
7k to 
75 St 


Kd 5 

Kell 

Kd 7 

d60 


Kell 
Kd 7 
g 36 
Kd 7 
e 36 
Kell 

76 tr d60 


yo . ^ 

y —^ y_position^ j find coordinate of first 
of first grid line horizontal grid line. 
Ay/2 1 position has heen fonnd 

when ovei*flow occurs. 


\ 

\56 

\57 

SU 

to 

tr 

Kd 7 

d60 

d52 

A y /2 J 

:6,d76) 

■>d60 

ca 

KdO 

1.00077 

^61 

st 

KclO 

X -^0 .77700 j 

[62 

ca 

Kell 

y “ coordinate / 

63 

ds 

KclO 

of grid line. > 
x-^ 0.77700 1 

61 ^ 

ca 

KclO 

x-» 0.77700) 1 

65 

ad 

Kd 3 

Ag J 

\ 66 

/to 

d70 

\ 67 / tr 

d6l 







display horizontal 
grid line. 


Ay/2 

Ay/2 

y-fAy 


find coordinate of next 
horizontal grid line. 

all horizontal lines have 

heen displayed when over¬ 
flow occursj at that time, 
transfer to leave. 


if horizontal grid incomplete, return to 
display next line. 


KdO 1.00077 
1 

2 0.—00 

3 0.—00 

4 

5 

6 0.-—0 

7 0.^—0 


H|no. of times to display axes. 

A.: distance between points in an axis. 

Ag S distance between points in a grid line. 


4 

Ay /2 / 


yo 


coordinates of origin< 
half spacings between grid lines. 





22-SA-56 


Page 19 


(e) SPECIAL SQUARE ROOT 


Input ; Enter at eO with the output of the SQUARE ROOT subroutine 
(SRL of memorandiaa 6 m- 3^97) AC, having previously- 
stored the scale factor (n) of x in KeO. (See Notes.) 

Output ; -h'fiT .2'“^ in AC, when x is stored as x .2”^. (See Notes) 

Other Subroutines Needed ; SRI#7f assiomed to be stored in field k. 


Duration; 159 microseconds, including the time spent after re-entering 
SRL # 7 . 

Accuracy; Maximum error is O.OOO 9 O (decimal), including 

that possible in the original ou-cput of SRL ^ —- this 
error requiring a scale factor of 2“ii, the same as the output. 


Storage; h program registers, 1 storage, plus l4 (decimal) regis-ters 
of SRL # 7 . 


Notes; 1. If n^O, the result obtained from use of the SCALE FACTOR 

subroutine (SRL ^ of memorandum 6 m- 3^9T) possesses a 
different scale factor than that of x. The subroutine 
des cribed here takes the results of SRL if^ and produces 
■ f x - with the original scale factor of x. For example; 

may be desired to app ear as 0 . 00002 ; whereas SRLl #7 
gives 0.00552, which is 2 " 13 , 

2, As used here, n “-^-15 indicates a scale factor of 2"15^ 
or, in other words, that the binary point is. considered to 

be at the right-hand end of the register 

3 . Each relative address denoted with a g in the following 
program refers to the decimal address as given in SRL ^7 
of memorandum 6 m-3497- 


Program; 


eO 

rf 

g 28 

prepare to return to main program. 

1 

st 

Kgl 

'fx', as obtained from SRL 

2 

ca 

KeO 

n 

3 

tr 

gl7 

transfer to last part of SRL ^0 


KeO -- n; original scale factor of x. 


J.%SKAVITCH 
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